Cosmography with the Einstein Telescope 
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Einstein Telescope (ET) is a 3rd generation gravitational-wave (GW) detector that is currently 
undergoing a design study. ET can detect millions of compact binary mergers up to redshifts 2-8. A 
small fraction of mergers might be observed in coincidence as gamma-ray bursts, helping to measure 
both the luminosity distance and red-shift to the source. By fitting these measured values to a 
cosmological model, it should be possible to accurately infer the dark energy equation-of-state, dark 
matter and dark energy density parameters. ET could, therefore, herald a new era in cosmology. 
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PACS numbers: 04.30.Db, 04.25.Nx, 04.80.Nn, 95.55.Ym 



The goal of modern cosmology is to measure the ge- 
ometrical and dynamical properties of the Universe by 
projecting the observed parameters onto a cosmological 
model. The Universe has a lot of structure on small 
scales, but on a scale of about 100 Mpc the distribu- 
tion of both baryonic (inferred from the electromagnetic 
radiation they emit) and dark matter (inferred from large 
scale streaming motion of galaxies) components is quite 
smooth. It is, therefore, quite natural to assume that 
the Universe is homogeneous and isotropic while describ- 
ing its large-scale properties. In such a model, the scale 
factor a(t), which essentially gives the proper distance 
between comoving coordinates, and curvature of spatial 
sections fc, are the only quantities that are needed to fully 
characterize the properties of the Universe. The metric 
of a smooth homogeneous and isotropic spacetime is 



da'' 



1 - fcCT2 



'■OdLp') 



where t is the cosmic time-coordinate, (ct, 9, Lp) are the 
comoving spatial coordinates, and /c is a parameter de- 
scribing the curvature of the t = const, spatial slices. 
k = 0, ±1, for flat, positively and negatively curved 
slices, respectively. The evolution of a{t) depends on the 
parameter fc, as well as the "matter" content of the Uni- 
verse. The latter could consist of radiation, baryons, dark 
matter (DM), dark energy (DE), and everything else that 
contributes to the energy-momentum tensor. 

The Friedman equation, which is one of two Einstein 
equations describing the dynamics of an isotropic and ho- 
mogeneous Universe, relates the cosmic scale factor a{t) 
to the energy content of the Universe through 



H{t) = Ho 



1/2 



(1.1) 



where H{t) = d{t)/a{t) is the Hubble parameter {Hq = 
H{tp) being its value at the present epoch tp), while 
Cluit) and ClA{t) are the (dimensionless) energy densi- 
ties of the DM and DE, respectively. The above equa- 
tion has to be supplemented with the equation-of-state 
of DM, assumed to be pressure-less fluid p = [tluit) = 
^m{^ + z)^, where J^m = ^M{tp)] and of DE, assumed to 



be of the form p ~ wpt^ [i^A(i) = ^tSX + zf''^^^^^ ^ where 
r^A = f^A(ip)], with w = —1 corresponding to a cosmo- 
logical constant. The goal of cosmography is to measure 
(i/o, f^M, f^A, w, k, . . .), which essentially determine the 
large-scale geometry and dynamics of the Universe. In 
the rest of this paper we shall assume that the spatial 
slices are flat (i.e., k = 0). 

Astronomers use "standard candles" to measure the 
geometry of the Universe and the various cosmological 
parameters. A standard candle is a source whose in- 
trinsic luminosity L can be inferred from the observed 
properties (such as the spectral content, time- variability 
of the flux of radiation, etc.). Since the observations 
also measure the apparent luminosity F, one can de- 
duce the luminosity distance Dl to a standard candle 
from _Dl = \/L/{4:TtF). In addition, if the red-shift z to 
the source is known then by observing a population of 
such sources it will be possible to measure the various 
cosmological parameters since the luminosity distance is 
related, when fc = 0, to the red-shift via 



A, = 



c{l + z) 
Ho 



dz' 

[r!M(l + ^')'^ + ^^A(l + ^')^^'+'"^ 



1/2- 

(1-2) 

There is no unique standard candle in astronomy that 
works on all distance scales. An astronomer, therefore, 
builds the distance scale by using several steps, each of 
which works over a limited range of the distance. For in- 
stance, the method of parallax can determine distances 
to a few kpc, Cepheid variables up to 10 Mpc, the TuUy- 
Fisher relation works for several tens of Mpc, the Dn-cr 
relation up to hundreds of Mpc and Type la supernovae 
up to red-shifts of a few [T]. This way of building the 
distance scale has been referred to as the cosmic distance 
ladder. For cosmography, a proper calibration of the dis- 
tance to high red-shift galaxies is based on the mutual 
agreement between different rungs of this ladder. It is 
critical that each of the rungs is calibrated with as little 
an error as possible. 

Cosmologists have long sought for standard candles 
that can work on large distance scales without being de- 
pendent on the lower rungs of cosmic distance ladder. In 
1986, one of us pointed out |2] that gravitational astron- 
omy can provide such a candle, or, more appropriately. 
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FIG. 1: The left panel shows the range of the Einstein Telescope for inspiral signals from binaries as a function of the intrinsic 
(red solid line) and observed (blue dashed line) total mass. We assume that a source is visible if it produces an SNR of at least 
8 in ET. The right panel shows a realization of the source catalogue showing the measured luminosity distance (inferred from 
GW observation of neutron star-black hole mergers) versus their red-shift (obtained by optical identification of the source). 
This catalogue is then fitted to a cosmological model. 



a standard siren, in the form of a chirping signal (i.e., a 
signal whose frequency increases as a function of time) 
from the coalescence of compact stars (i.e., neutron stars 
and black holes) in a binary. The basic reason for this is 
that the gravitational-wave (GW) amplitude depends on 
the ratio of a certain combination of the binary masses 
and the luminosity distance. For chirping signals GW 
observations can measure both the amplitude of the sig- 
nal and the masses very accurately and hence infer the 
luminosity distance. 

Let us first recall in some detail how we might mea- 
sure the luminosity distance. We will first assume that 
the source is located close-by, i.e., its redshift z ^ 1, 
although we will later relax this condition. 

Gravitational waves are described by a second rank 
tensor h^p, which, in a suitable coordinate system and 
gauge, has only two independent components /i+ and , 
hxx = ~hyy — h^, hxy = kyx = /ix, aU othcr compo- 
nents being zero. A detector measures only a certain 
linear combination of the two components, called the re- 
sponse h(t) given by 

hit) = F+{e, ^, ij)h+it) + (e, if, ij)h^ {t), (1.3) 

where F+ and F^ are the detector antenna pattern func- 
tions, tjj is the polarization angle, and (6', t^) are angles 
describing the location of the source on the sky. The an- 
gles are all assumed to be constant for a transient source 
but time-dependent for sources that last long enough so 
that the Doppler modulation of the signal due to the 
relative motion of the source and detector cannot be ne- 
glected. For a coalescing binary consisting of two stars 
of masses mi and m2 (total mass M = mi + m2 and 
symmetric mass ratio v = mim2/M^) and located at a 
distance Dl, the GW amplitudes are given by 

h+{t) = 2iyM^/^ D^\l + cos^{L))uj{t^tQf/^ 

X cos[2<t>{t - to; M, i^) + ^a], 

hy,{t) = AvM^'^ D:[^cos{i)uj{t~toY'^ 

X sin[2$(t-to; M, v)+'^'ol (1-4) 

where b is the angle of inclination of the binary's orbital 
angular momentum with the line-of-sight, <jj{t) is the an- 
gular velocity of the equivalent one-body system around 



the binary's centre-of-mass and $(t; M, v) is the corre- 
sponding orbital phase. Parameters to and 4>g are con- 
stants giving the epoch of merger and the orbital phase 
of the binary at that epoch, respectively. 

The above expressions for /i+ and h-^ are the dominant 
terms in what is essentially a FN perturbative series. We 
have written down the expressions for a system consisting 
of non-spinning components on a quasi-circular orbit. In 
reality, we cannot assume either to be true. Eccentric- 
ity might be negligible only in the case of stellar mass 
binaries expected to be observed by ground-based detec- 
tors, but both eccentricity and spins could be non-zero in 
the case of merger of supermassive black holes expected 
to be observed by the Laser Interferometer Space An- 
tenna (LISA) [3^. The argument below holds good to 
whatever order the amplitudes are written down and for 
non-spinning objects on an eccentric orbit. 

Substituting the expressions given in Eq. (1.4 1 for 
and ft.x in Eq. (1.3 1, we get 



h{t) = 



,2/3 



cos[2$(< - to\ M, v) + «>o],(1.5) 



eff 



[F2 (1 + COs2(t))2 + 4F2 C0S2(/,)] 
2i^x COs(t) 



<i>n -I- arctan 



i^+(l-HC0s2(t)) 



172' (1-6) 
(1.7) 



Here Doff is the effective distance to the binary, which 
is a combination of the true luminosity distance and the 
antenna pattern functions. Note that D^s ^ ^l- In the 
case of non-spinning binaries on a quasi-circular orbit, 
therefore, the signal is characterized by nine parameters 
in all, {M,i',to,^o,0,'P,i',>',Di,)- 

Since the phase ^{t) of the signal is known to a high 
order in FN theory, one employs matched filtering to ex- 
tract the signal and in the process measures the two mass 
parameters (M, ly) (parameters that completely deter- 
mine the phase evolution) and the two fiducial parame- 
ters {to, <I>o)- In general, the response of a single inter- 
ferometer will not be sufficient to disentangle the lumi- 
nosity distance from the angular parameters. However, 
EM identification (i.e., electromagnetic, especially opti- 
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FIG. 2: The plot on the left shows the distribution of errors in Qm, and w, obtained by fitting 5,190 realizations of a 
catalogue of BNS merger events to a cosmological model of the type given in Eq. (1.21, with three free parameters. The 
fractional 1-a width of the distributions (jQ^/riM, crn^/SlA, and cr„/|w|, are 18%, 4.2% and 18% (with weak lensing errors in 
-Dl, left panels) and 14%, 3.5% and 15% (if weak lensing errors can be corrected, right panels). The plot on the right is the 
same, but assuming that Qa is known to be Ha ~ 0.73, and fitting the "data" to the model with two free parameters. The 
fractional 1-a widths in the distribution crn^/nM and and auj/|w|, are 9.4% and 7.6% (with weak lensing errors in Dl, left 
panels) and 8.1% and 6.6% (if weak lensing errors can be corrected, right panels). 



cal, identification) of the source will determine the direc- 
tion to the source, still leaving three unknown parameters 
L, Dl)- If the signal is a transient, as would be the 
case in ground-based detectors, a network of three inter- 
ferometers will be required to measure all the unknown 
parameters and extract the luminosity distance. 

Although the inspiral signal from a compact binary is 
a standard siren, there is no way of inferring from it the 
red-shift to a source. The mappings M (1 + z)M, 
uj ujI{1 + z), and Dl ^ (1 -I- z)Dl, in Eq. ( [Til , leave 
the signal invariant. Note that a source of total mass 
M at a red-shift z will simply appear to an observer to 
be a binary of total mass (1 -I- z)M. One must optically 
identify the host galaxy to measure its red-shift. Thus, 
there is synergy in GW and EM observations which can 
make precision cosmography possible, without the need 
to build a cosmic distance ladder. 

Over the next two decades GW interferometric de- 
tectors will provide a new tool for cosmology. Ad- 
vanced ground-based interferometers, operating around 
2015, are expected to detect ^ 40 binary neutron star 
mergers each year from within about 300 Mpc. Redshift 
could be measured to a (small) number of events asso- 
ciated with GRBs, thereby allowing an accurate deter- 
mination of the Hubble constant [U [S] . Observation by 
LISA of extreme mass ratio inspirals could measure the 
Hubble constant pretty accurately j6j. LISA will also 
observe binary super-massive black hole mergers with 
SNRs ^ few X 1 , 000 enabling the measurement of the 
DE equation-of-state to within several percents [7, 8J. 

In the rest of this paper we will discuss how well it 
might be possible to constrain cosmological parameters 
by GW observations of the inspiral signal of compact bi- 
naries by the Einstein Telescope (ET) — a third genera- 
tion GW interferometer that is currently under a design 
study 0. ET is envisaged to be ten times more sensi- 
tive than the advanced ground-based detectors, covering 
a frequency range of 1-10"* Hz, posing new challenges in 



mitigating gravity gradient, thermal and quantum noise. 

The sky-position averaged distance up to which ET 
might detect inspiral signals from coalescing binaries with 
an SNR of 8 is shown in Fig. [l] We plot the range both as 
a function of the intrinsic (red solid lines) and observed 
(blue dashed lines) total mass. A binary comprising two 
1.4 M0-neutron stars (BNS) can be observed from a red- 
shift of z ~ 2, and that comprising a 1.4M0-neutron star 
and a 10 Mg-black hole (NS-BH) from z ~ 4. 

The expected rate of coalescences per year within the 
horizon of ET is - several x 10^ for BNS and NS-BH. 
Such a large population of events to which luminosity dis- 
tances are known pretty accurately, would be very useful 
for measuring cosmological parameters. If, as suspected, 
BNS and NS-BH are progenitors of short-hard gamma- 
ray bursts (GRBs) [lOj, then it might be possible to make 
a coincident detection of a significant subset of the events 
in GW and EM windows and obtain both the luminosity 
distance to and red-shift of the source. 

Since GRBs are beUeved to be beamed with beaming 
angles of order 40° , we assume that only a small fraction 
(~ 10"'^) of binary coalescences will have GRB or other 
EM afterglows that will help us to locate the source on 
the sky and measure its red-shift. Eventually, we will 
be limited by the number of short-hard GRBs observed 
by detectors that might be operating at the time. As a 
conservative estimate, we assume that about 1,000 BNS 
mergers will have EM counterparts over a three-year pe- 
riod. For definiteness we consider only BNS mergers and 
take these to have component masses of (1.4, 1.4)M0. 

How well would we measure cosmological parameters 
with a catalogue of such sources? To answer this ques- 
tion we simulated 5,190 realizations of the catalogue con- 
taining 1,000 BNS coalescences with known red-shift and 
sky location, but the luminosity distance subject to sta- 
tistical errors from GW observation and weak lensing. 
One such realization is shown in Fig. [T| (right panel). 
We assumed that the sources were all in the red-shift 
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range < z < 2, distributed uniformly (i.e., with con- 
stant comoving number density) throughout this red- 
shift range. The luminosity distance to the source was 
computed by assuming an FRW cosmological model with 
Hq = 70kms-iMpc-i, = 0.27, f^A = 0.73, and 
w = — 1, but the measured distance was drawn from a 
Gaussian distribution whose width ct^Jl was determined 
by the quadrature sum of the errors due to weak lens- 
ing and GW observation. Weak lensing error in 
was assumed to be 5% at z = 1 and linearly extrapo- 
lated to other red-shifts. GW observational error was 
estimated from the covariance matrix Ckm of the five- 
dimensional parameter space of the unknown signal pa- 
rameters pk = (M, to, ^"0, £'l): 



A 



A. 



{hk, hj, 



dh 
dpk 



(1.8) 



Here the angular brackets denote the scalar product, 
which, for any two functions ait) and bit), is defined as 



(a, b)^m f 
Jo 



df 

Shif) 



A{f)B*{f) 



(1.9) 



where A and B are the Fourier transforms of the func- 
tions a(t) and h{t), respectively, and Sh{f) is the ET 
noise power spectral density. Note that since GRBs are 
expected to be strongly beamed, we did not take the an- 
gles {i,4') associated with the unit normal to the plane 
of the inspiral as unknown variables. This assumption 
is justified: even if the opening angle of a GRB beam 
is as large as 40°, the unit normal to the plane of the 
inspiral would still be confined to only 3% of the area of 
a unit sphere. Averaging errors over (t, ip) with the con- 
straint i < 20° would then be little different from taking 
i = 0°. We did, however, average the errors over the sky 
position angles (6*, </>). We then fitted each realization 
of the source catalogue to the cosmological model given 
in Eq. (1.2|, using the Levenberg-Marquardt algorithm 



[TTl [T^ . in order to find a set of best fit parameters. It 
turns out that a catalogue of 1,000 sources is not quite 
enough for an accurate determination of all the parame- 
ters. However, assuming that Hq is known accurately, the 



algorithm gave the best fit parameters in {^Im, ^a, w) for 
each of the 5,190 realizations. 

The distribution V of the parameters obtained in this 
way are shown in Fig. [2] where the vertical line is at the 
true value of the relevant parameter. The relative 1-a 
errors in Hm and w, are 4.2%, 18% and 18% (with 
weak lensing) and 3.5%, 14% and 15% (with weak lensing 
errors corrected). Although Vi^w) is quite symmetric, 
V{flM) and V{flA_) are both skewed and their mean values 
are slightly off the true values. The medians, however, 
are coincident with the true values. 

In addition to Hq if f^A is also known (or, equivalently, 
if ilivi + = 1), then one can estimate the pair (f2M, w) 
more accurately, with the distributions as shown in Fig. [2] 
with greatly reduced skewness, and 1 -fj errors in and 
w, of 9.4% and 7.6% (with weak lensing) and 8.1% and 
6.6% (with lensing errors corrected). Finally, if w is the 
only parameter unknown, it can be measured to an even 
greater accuracy with 1-a errors of 1.4% (with weak lens- 
ing) and 1.1% (with lensing errors corrected). 

The results of our simulation are quite encouraging but 
further work is needed to confirm the usefulness of GW 
standard sirens in precision cosmology. Let us mention 
some that are currently being pursued. Spins of com- 
ponent stars can be legitimately neglected in the case 
of neutron stars (and hence in BNS) but not for black 
holes. The modulation in the signal caused by the spin 
of the black hole can improve parameter accuracies. We 
assumed, for simpHcity, that all our sources are BNS sys- 
tems with masses (1.4, 1.4)M0. In reality, the catalogue 
will consist of a range of NS and BH masses. A more re- 
alistic Monte Carlo simulation would draw binaries from 
the expected population rather than the same system, 
some of which (e.g. more massive systems) would lead 
to better, but others to worsened, parameter accuracies. 
The signal contains additional features, such as other har- 
monics of the orbital frequency than the second harmonic 
considered in this work and the merger and ringdown sig- 
nals. These are important for heavier systems and could 
potentially reduce the errors. These factors are currently 
being taken into account to get a more reliable estimation 
of the usefulness of ET in precision cosmography. 
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